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Abstract 

We describe some ways how higher order corrections can reveal themselves if integrated 
over the infrared region. We show that in different renormalization group (RG) schemes 
Q\ ', and for some observables one has no factorial divergences. We argue that for treating 

things in the infrared region it is preferable to start with a RG scheme without the 
■ infrared Landau pole in the running coupling constant. The uncertainties for the r lepton 

|"t j ' width resulting from accounting for higher order corrections are discussed. 
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With new high order corrections of perturbation theory hardly available anymore in 
cases like e + e~ annihilation or r lepton width |]J it is tempting to speculate on the 
general structure of series within perturbation theory (PT) |3|, |J. Some attention 



has been recently paid to possible factorial divergences in PT series for observables that 

O ■ 

include integration over an infrared region in momentum space [gj, || [7|, |S|. At the level 
of diagrams within PT it happens to any observable and the factorial divergence due 
to the simple bubble chain diagrams is predicted. In fact this prediction is not well 
Q_i! justified in QCD because of the choice of a particular subset of diagrams and even of 

the special contributions (nonabelianization ||). Next terms of PT expansion can cancel 
these divergences because such terms become large at small momenta and can not be 
^ ■ treated as corrections. In QED there is a formal parameter for ordering diagrams within 

the 1/Nf expansion and the statement about the factorial growth of coefficients can be 
confirmed by direct computation [|IIJ but it is known also that QED within the 1/Nf 
expansion has practically no sense at all and using it as a guide for general structure of 
QCD is not well grounded. 

In this paper we consider some observables that are represented by integrals over the 
infrared region and give several ways to define them using the freedom of choice of the 
RG scheme. The main conclusion we draw is that the results of integration can easily be 
made well defined without any nonperturbative (in a strictly defined sense) contributions. 
These results are ambiguous to the same degree as any ordinary PT series, numerically it 
can be important because in the infrared region the coupling constant becomes large in 
most of "natural" RG schemes. However it can be made small as well by some particular 
choice of extrapolation to low momenta. 

The paper is organized as follows. First we define our main object as an integral 
of some PT expansion over infrared region and show that the corresponding expression 
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can contain no factorials depending on the integrand. Even with the simplest version of 
RG analysis including only the leading (3 function and summation of leading log's, the 
Landau ghost contributions to different orders in a s can cancel each other. This leads 
us to conclusion that infrared renormalons or factorial growth of coefficients due to some 
infrared integration has no invariant meaning. 

Next we present a number of examples where integrals are well defined and discuss 
their properties as functions of the boundary value of the coupling a s (s). It is shown that 
integrals can be analytic or nonanalytic in a s (s) in the vicinity of a s (s) = depending on 
a (3 function used. So, extrapolation of PT evolution to low energy domain can be done 
in different ways and the results vary strongly. 

These two parts accomplish our theoretical consideration of the problem with conclu- 
sion that nothing definite can be said about the structure of higher order terms on the 
basis of some subset of diagrams. 

Then we go to practical applications and estimate uncertainties of theoretical predic- 
tions for the t lepton width when different approaches for defining the integration over 
infrared region are used. This problem has been widely discussed in the literature (as 
recent references see, e.g. [[II], [H| so we limit ourselves to qualitatively different 
versions of changing RG schemes only. The basic (or naive) scheme is one where only 
self-consistent number of terms of the expansion in the coupling constant is kept. Note 
that adding all terms proportional to RG log's can produce any answer - just the re- 
flection of the Landau singularity. Then we use the K scheme JEIL [T5|, [TB| where there 



are no corrections to R(s) and the corresponding (3 function in three-loop approximation 
acquires an infrared fixed point so that the coupling constant can be extrapolated to the 
origin without any singularity. 

Finally we propose a set of schemes that regularize the infrared behavior of the coupling 
constant in general and allow one to use any reference scheme for high energy domain. All 
these schemes are legal and perturbatively equivalent at high energies. The uncertainties 
that come from low energy region are quite essential as our study shows. 

Thus the main object of our interest is the expression of the form 

F(s) = - f p(t)u(t)dt (1) 



s Jo 

where p(t) is given by a PT expansion 

p(t) = a s (t) + ... 

while u(t) is a kinematical factor smooth in the integration region. Expression ([I]) enters 
the finite energy sum rules analysis flT7| ]. We do not dwell upon the peculiarities connected 
with the analytic continuation [pl|, |l|J and consider eq. ([!]) as a formal expression. Note 
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that the formulation with the help of integration along the circle in the complex plane 
does not help much by itself because RG improved expressions for PT approximations 
of physical quantities have wrong analytical properties so the spurious singularities give 
contributions to physical results p0| . 

We restrict ourselves to the case of u(t) — 1 and compare 



1 r s 

F(s) = - p(t)dt and p(s) = a s (s) + 
s Jo 



For technical convenience we introduce a new notation for the coupling constant a(s) = 
(3 a s (s) so that asymptotically a(s) — > (/ra(s/A 2 )) -1 when s — ► oo. A RG equation for 
our rescaled coupling constant a(s) reads 

s-^a(s) = /3(a) = -a 2 (l + ca 2 + da 3 + . . .) (2) 

with rescaled values of coefficients c = 64/81, C\ = 3863/4374 for rif = 3 in the MS 



renormalization scheme |21 



Taking p(t) = a(t) we find the factorial growth of coefficients of a(s) expansion of F(s) 
that reflects the presence of infrared Landau ghost (or divergence of the formal series for 
the integrand) 

F(s) = - / a(t)dt = a{s) V a n {s)n\, a{s) > (3) 
where in the leading approximation 

a w- i-<w) ^ a(s) I a " w '"" (s/i) - (4) 

Last equality in eq. (§) is formal and valid only within the convergence circle of the power 
series \a(s)ln(s/t)\ < 1, the use of this equality outside the convergence circle leads to 
all problems of factorial growth. Factorials can be canceled by higher order terms in 
the integrand p(t). Indeed, for p{t) = a(t) + ka 2 (t) we find for a simplified (3 function 
f3(a) = —a 2 with c = that corresponds to n/ ~ 8 

-/ (a(t) + ka 2 (t))dt = a(s) + a(s)(l + k)J2a n (s)n\. (5) 

S ^° ' n=l 

If k — — 1 there are no factorials in the expression of the form and only first term sur- 
vives. Thus higher order corrections to the integrand can drastically change the situation 
with factorials. It is natural because factorials come from the region where higher order 
corrections are not small and comparable (even much larger) than leading contributions. 
In general, we find for integrals of the invariant charge given by eq. ([D 
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and for p(t) = a(t) J2 p =oPp aP ( t ) 



1 r s N o 00 

- / p{t)dt = a{s)Y, P ^Y, a n (s)n 

b J0 p =0 n=p 



00 N N N-l 

= a(s) £ a»(*)n! £ ^ + «00 E ^ E « B («)n!. 

n=7V p=0 p=0 n=p 

The last addendum is a finite sum (an analog of a(s) in eq. (HD) while the coefficient in front 
of the infinite sum with factorial terms is a linear combination made from the coefficients 
p p of the observable p(t). Inclusion of nonvanishing coefficients of higher order terms of 
the /3 function is straightforward and only leads to extended combinatorics leaving the 
statement about the factorial growth unchanged if a (3 function keeps to be negatively 
defined for all a. Even if the factorial growth persists (not complete cancelation between 
contributions from different orders) its coefficient (the strength) is determined by a linear 
combination of coefficients p n coming from all orders. It is a consequence of the fact 
that at low energy the formal ordering with a(s) is invalid and all orders contribute a 
comparable amount to the coefficient of factorial. In other words, factorials appear due 
to RG log's but in higher orders coefficients of these log's contain not only the f3 function 
coefficients but lower order terms of observables themselves. So the behavior at large 
orders of PT for the integral can not be fixed without knowing the whole series for the 
p(t) that is hopeless. The finite series for a (3 function in schemes like 't Hooft's one is 
illusion - it is still a nonperturbative statement about vanishing coefficients in all orders. 

Now we go over to examples without the Landau ghost in the infrared region that are 
free of any problems with integration or factorial growth. The purpose can be achieved 
with different means. In higher orders a j3 function can have a fixed point in some scheme. 
This happens for e + e~ annihilation in third order in K scheme without corrections. Or 
one can write down some new f3 functions that are defined at all values of coupling and 
determine its evolution for all momenta in a smooth way. 

We start with the second option and consider the model (3 function EO 



,2 



/3(a) = -——, k>0 (6) 

1 + KCL Z 



while f3 as (a) = —a 2 + . . . Then the RG equation (^j has a solution 



-In-^ + ^ln^+AK 
<z) = (7) 

and the unphysical pole at z = A 2 of the asymptotic solution a as (z) = (ln(z / A 2 ))^ 1 dis- 
appears. Thus, the particular way of summing an infinite number of specific perturbative 
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terms for the (3 function can cure the Landau pole problem. Let us stress that there are 
no nonperturbative terms added but the freedom of choosing a renormalization scheme 
for an infinite series was used instead. This is however beyond the formal framework 
of perturbation theory where only finite order polynomials in the coupling constant are 
allowed as expressions for any quantity. 

Formula @ can be considered either as a pure PT result in some particular RG after 
an infinite resummation or as a sort of Pade approximation of some real (3 function that 
might include nonperturbative terms as well. The only important point for us here is 
that the running coupling obeying the RG equation with such a (3 function has a smooth 
continuation to the infrared region. 

The effective charge given by eq. (0) has a correct asymptotic behavior at z — > oo 
and no singularity in the whole complex plane with a cut along the positive semiaxis 
if one considers it as a running coupling constant in Euclidean space. These properties 
make it a good expansion parameter for physical observables that obey the dispersion 
relation because it has a proper analytic behavior unlike the asymptotic charge with the 
Landau pole. Note that this scheme does not solve the problem of strong coupling and 
even fails to bypass this problem when the expansion for an observable is alternate and 
higher order terms can make it negative at small enough s that can contradict spectrality 
for observables of the type of cross-sections. Because the expansion parameter becomes 
large in infrared region the polynomial approximation is invalid in this domain. Here 
we encounter a particular case of the general situation that the expansion in unphysical 
parameter a s is incorrect and the proper way of action is to expand one physical quantity 
through another. Such formulation of PT is physically more justified. The use of K 
scheme is a particular example of this approach. 

Now we turn to the consideration of the structure of expressions of the form (|l|) with 
different extrapolation of the running coupling to infrared region. In general, the RG 
equation for 



and the solution of this differential equation can be found with quadratures because it is 
a linear differential equation of the first order. At simple (3 functions like (3(a) = —a 2 the 
solution of the RG equation (an integral curve) fixed by the boundary condition a(s) — > 
at s — ► oo can not be continued to the point s = 0. The integral curve goes only until 
s = A 2 where there is a vertical asymptote. New information is necessary for continuation 




has the form 




(8) 
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of the solution to the origin. Consider now an example with 

-a 2 

3(a) = , a > 0. (9) 

In this case the interval (0, oo) in s is mapped uniquely to the interval (0, oo) in a and 
the solution can be uniquely continued to the origin. Still a special consideration of 
small s region is necessary to determine the correct contribution of the solution of the 
homogeneous equation. We find 

F(a) =a + a 2 + AF(a) 
where AF(a) is a solution of the homogeneous equation 

F(a) + 13(a) F' (a) = (10) 

AF(a) = F a 2 exp(--). 

a 

Coefficient F Q is determined from the infrared region to be F = — 1. The result is 

F = a + a 2 — a 2 exp( — ). 

a 

This result can be obtained by explicit integration as well. Namely, the RG equation for 
the effective charge of eq. (0) is given by 

^ 1 



ln(s/A 2 ) = - - Una 



a 



and F(s) can be rewritten as (ao = a(s)) 



F(s) = F(ao) = e- { ^- 2lnao) [°° ae^~ 2lna \- + 2^)da 



-(-L-2ina ) f 1/a ° 



r 1 / a o i 
/ eZ(£ + 2)d£ = a + a 2 (l-e~^). 
Jo 

The last term gives the "condensate" contribution to F(s) at a — >• (s — > oo). Up to 
logarithmic corrections 



F cond (s) ~ 



A2 

s 

At small s the effective charge runs as a(s) ~ A/y^ and 

2A 

F ( a r . 

Note that the 3 function @ does not regularize the whole physical quantity: higher order 
terms in the expansion (like a 2 ) can be nonintegrable. In this particular example it is the 
case. 
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This simple example shows that for a reasonable definition of the evolution in the 
infrared region factorials are absent in the expansion. The last term cannot be detected 
if the equation (H) is integrated by series near a = 0. Indeed, for F(a) = J2k=i fka k one 
has for eq. (§) with (3 function (]|) 

fx = l, f 2 = l, f k - (k - 3)A_! = 0, k>2. 

Thus, fk = 0, k > 2. Note that the last equation would lead to factorial growth but 
instead all higher order coefficients vanish. 

To study higher orders of PT the Borel transformation is often used. The Borel 
analysis for this example goes as follows 

F(a)= /°°e-ls(0de (H) 

and the Borel image B(x) is 

B(x) = l + 6(x-i)+x6(l-x). (12) 

So the PT series with all coefficients known still does not allow to restore the exact 
answer through Borel summation. A naive Borel image for polynomials from definition 
(PI) behaves as B(x) = Z)fc=i fk x /(& — 1)' that is correct at small x and completely 
wrong at large x and has nothing to do with the exact one (0). Being rewritten as 

B(x) = l + x + {l-x)6(x-l) (13) 

the result shows that Borel image is singular (nondifferentiable) at x — 1 though this is 
not a pole singularity. 

The corresponding example with f3 function (^) can be written as follows 

F(s) = a + -a 3 = - / (a - a 2 + -a 3 )ds. 
3 s Jo 3 

The RG equation of the type (§) can be integrated by series but the homogeneous equation 
(PH|) has no solution satisfying the boundary conditions at small s (constraint on the 
growth). To understand this feature we note that at small s the running coupling behaves 
like a(s) ~ K^ 1 ln(A/ s) and after integration the logarithmic behavior survives. As for the 
solution of the homogeneous equation it contains a factor 

a 

that is not logarithmic at small s (a — > oo) and can not appear in the resulting expression. 
In contrast to the previous example all powers of the effective coupling are integrable here. 
The rate of increase is however quite large for higher order terms. 
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Now we define a set of schemes that provide a smooth extrapolation to the infrared 
region that are a direct generalization of the model (||). The k scheme is determined by 
the (3 function 

«•) = T3§b (14) 

where (3(a) is a (3 function in a reference scheme, MS for instance. The form flTj]) is 
chosen because of practical convenience only - it requires no more work as corresponding 
reference scheme and eliminates the Landau pole in infrared region. 

The (3 function given by eq. (|TJ]) is bounded at large a and eq. (Q) has a solution for 
a(z) that is defined on the whole positive semiaxis and is free of the Landau pole. The 
absence of singularities (Landau ghost) allows one to use the evolution of coupling till 
the origin and destroys the mathematical part of the reasoning of ref. [|5| about Borel 
nonsummability and also Parisi's conjecture |E| looks wrong in this particular scheme. 

The solution of the type (|7|) can be rewritten without using A in a more traditional 



form through some intermediate energy scale ji The parameter A has no special 



meaning anymore as a position of the pole of the running coupling. Thus, in the pure 
PT framework one can eliminate the Landau pole of the coupling using the freedom of 
the renormalization scheme choice. The smooth coupling however does not require any 
nonperturbative effects to regularize the integrals in which it appears. 

The solution to the RG equation for the invariant charge in the k scheme is simple 
because it is closely related to the MS running coupling 

n+l 

Inis/K 2 ) = $(a) - k- (15) 

w ' w n+ 1 

where 

f doc 



(3(x) 



and the normalization is chosen in the form 



with f3(2)(a) = — a 2 (l + ca). This gives the standard definition of the parameter A to be 
A^jg. Two charges are connected through 

■ • K a n+3 + o(a n+3 ). 



n + l 



The exact connection can be found from eq. (|T5|). So, for n > 1 they coincide for all 
observables because in practice there are no calculation beyond the third order. 

Still if the corresponding expansion of a physical observable is alternative, higher order 
corrections can make the physical quantity negative in some part of the integration region. 
In this sense the approach does not give the full regularization. 
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Concluding the theoretical part of the paper we state that there are different ways 
to bypass the problem of factorial growth of coefficients for integrals from physical ob- 
servables. The simple change of RG scheme within PT in high energy region can allow 
a smooth extrapolation of the effective charge into low energy domain that makes all in- 
tegrals well defined. The change of scheme leads to numerical uncertainties in prediction 
that is a common feature of any PT result. We now consider these uncertainties for predic- 
tions of the r lepton width or for the parameter a s (m 2 .) extracted from experimental data 
on this width. We stress that we limit ourselves only to "qualitatively" different schemes 
that describe different infrared behavior of an invariant charge. We do not discuss small 
variation of schemes within the qualitatively equivalent class. 

The expression for the r lepton width has the form ]23 



Rr= l'" T ^2(l-^ni + 2^)R(s) (17) 
mz mz mr 



T 



with R(s) given by 



R(s) = 3(1- 



7T 

Using the experimental value (as a recent reference, see ||11|| ) 

R e T xp = 3.645 ± 0.024 (18) 

and writing 

Rr = 3(1 + £r T ) 
Po 

we find 

r r = a(/i 2 ) + . . . and rf p = .4838 ± 0.018. 

First, for the sake of completeness and for purposes of normalization we compute the 
naive RG result that corresponds to keeping only three terms of expansion in a(/x 2 ) with 
explicit RG logarithms. Having rewritten R(s) as 

J2( a ) = 3(l + £r( a )), 

one gets for properly normalized QCD contribution (we do not discuss any corrections 
due to electroweak interactions and take this process only for demonstration of theoretical 
uncertainties stemming from the choice of RG scheme in QCD) 

r(s) = a(/i 2 ) + aVX^i + L) + a\^ 2 ){k 2 + L(c + 2k x ) + L 2 ) (19) 



with k\ = 0.7288, ki = —2.0314 ]23, 2"5(| , L = ln^-. Introducing normalized moments 

N 

,2 
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we get the well known result 

r T = 2(r - r 2 + r 3 /2) 

and 



(N + l) f 
Jo 

This gives 



x N {—InxYdx - 1 



(N + 1)p 



rK = a + a ^ 1 + _j +a 3( fe + _ + __j. (20) 

Collecting all together we get the result 

r r = a + (ki + 5i)a 2 + (k 2 + (c + 2fei)<Ji + 5 2 )a 3 = a + k[a 2 + k' 2 a 3 (21) 

where 5i = 19/12, S 2 = 265/72. The equation for determination a s (m^) reads 

r T = r e T xp = 0.4838 ± 0.018. (22) 

We have two options here: to solve equation (^) with the theoretical expression for r T 
given by eq. (refresultth) exactly or perturbatively. Exact answer is 

a = 0.253 ± 0.006 and a s (m 2 T ) = 0.353 ± 0.008. 

Perturbative method gives 

a(m 2 T ) = TV - k[r 2 T - (k' 2 - 2(k' 1 ) 2 )r 3 T + 0(4) 

= r T - (ki + 5^1 - (k 2 + (c + 2k x )8 x + 5 2 - 2(h + S^r* + 0(4) = 0.564 ± 0.046 
and 

4:71 

a s {m 2 T ) = —a(m 2 T ) = 0.787 ± 0.063. 
9 

Both results are purely perturbative and the errors reflect only the experimental error 
with the theoretical procedure being strictly fixed. We do not consider any options with 
choice of RG without introducing any qualitative change. We also do not discuss any 
keeping additional terms known due to RG (like a™ln n (n 2 / s)). 

Our last remark here is to formulate the result in a RG invariant way in the spirit of 
ref. [26]. The simplest (not unique) way is to express the integral through the value of 
integrand taken at the end point. If r(s) = a + k x a 2 + k 2 a 3 then 

- f r(t)dt = r(s) + r(s) 2 + (c + 2)r(s) 3 (23) 
s Jo 

and for the r lepton width we find 

19 /19 265\ 

r r = r(s) + 5 ir (s) 2 + (S lC + 5 2 )r(s) 3 = r(s) + j^r(s) 2 + i^—c + —j r(s) 3 . (24) 
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Eq. (]24D corresponds to the spirit of PT in that respect that it gives the direct relation 
between two measurable quantities. As for numerical values, the exact solution is r(m^) = 
0.270 while the perturbative one is r PT (m^.) = 0.123. For corresponding couplings we find 
a PT {m 2 T ) = 0.164 and a e s xact (m^) = 0.359. So, these two observables are poorly connected 
by PT relation that means that the operation of integration disturbs the PT expansion 
for r(s) strongly. From formula (|23| ) and integration rules for log's we deduce that large 
contributions of infrared region should be suppressed, i.e. only computation of high order 
moments are allowed as a PT operation. For three orders of coupling constant expansion 
and for the Nth order moment we require that the change of coefficients would be small 
enough 

M0 = r + ^+(^ + _L p ) r '. (25) 

This relation (]2~5D does not include any coefficients k t of the expansion of r in unphysical 
parameter and reflects only the influence of integration. Requiring this influence 
to be perturbative we get iV > 1 — 2 for particular values of r and c. The estimate 
depends on the number of RG log's kept for the expansion: the more log's kept in the 
expansion the larger the low limit for moments remaining perturbatively related to the 
integrand. This example again shows that the MS coupling constant is artificial (and 
auxiliary) quantity. And with integration over regions where PT series is invalid the 
integration measure can not be chosen arbitrary. Of course, numerically this particular 
computation can be made more appropriate by choosing some intermediate point s* for 
normalization instead of the end point and by using r(s*), s* < vr? T as an expansion 



parameter for the integral. This method would correspond to one of ref. [p7 , p8f but 
for fixed order of moment s* approaches to A 2 if higher orders of log's would be kept. 
Still, such a formulation of the calculation is more physical one because it connects two 
measurable quantities directly without intermediate things like The validity of PT 

expansion in this case means that two quantities commensurate and can be a measure 
for each other. We see that integration through infrared region violates this property 
unless the weight function is properly chosen. The choice depends also on the order of 
PT: at higher orders the suppression must be larger. So, the result of integration and the 
integrand are poorly connected by PT. If we had only higher moments (say, third one) 
for the width the situation would be different for this particular order of PT. 

In fact, the most successful method in this framework would be the Stevenson's ap- 
proach PS|, 3U| that requires an independent optimization for the integral as an indepen- 
dent PT quantity. Being extremely flexible it creates the best scheme for every process 
and every order of PT independently. We did not check however whether this method 
gives the convergent answer for the problem at hand. 



11 



Next we consider methods that allow to integrate after applying the RG that is beyond 
the scope of conventional PT. 

First, in three-loop approximation in the K scheme there is a fixed point for the running 
coupling and consequently there is no Landau pole problem. Namely, in the scheme where 
there are no corrections to r(s) the running coupling reads 

r(s) = clk = a + k\a 2 + k 2 a 3 

while the (3 function in this scheme is 

/3k(cik) = —o>k (1 + ca K + (ci - cki - k\ + k 2 )a 2 K + ...) = —Or (1 + ca K + c[a 2 K + . . .) 

with c[ = —2.2552. Fixed point is given by the equation 

1 + CCLx + C\d 2 K — 

with the solution ax = 0.8637, and (a s )K = (47r/9)ax = 1.206. The relation between 
coupling in K and MS schemes is 

ax = a + kia 2 + k 2 a 3 , a = ax — k\a K — (k 2 — 2k\)a K . 

Then the procedure is as follows. We compute integrals exactly an find 

a s (m 2 T ) K = 0.380 ±0.013. (26) 
Computation of MS coupling constant gives two answers again 

a s {m 2 T ) exact = 0.361 ±0.014 

a(m 2 T r l = a K (m 2 T ) - ha K (m 2 T ) - (k 2 - 2fc?)<&(rr£), 

a s (m 2 T ) po1 = ^a{m 2 T Y° l = 0.392 ± 0.017. 
9 

Note that the difference between them is larger than the uncertainty due to the experi- 
mental error and two intervals hardly overlap. 

Finally, we give the results in our k scheme analysis that is more general and universally 
applicable to integration in any basic scheme as well. This is a generalization of the fixed 
point approach. 

The spectral density in k scheme (with n = 0) up to third order is 

p(s) =a K + k x a\ + (k 2 + K)a 3 K , 
a K = a — kcl 3 , a = a K + Ka 3 K . 
12 



The expansion of (3 function reads 



/3{a K ) = -a 2 K (l + ca K + (c t - K)a 2 K + . . .) 

so introducing k is equivalent at high energies to change of c\. At low energies however 
they are different. 

Integration for moments can be easily rewritten in term of the charge itself 



m r ds , s 



r N = (N + l) ^(—r P (a) 
Jo m; m A T 



= (N + 1) P exp(N + 1)($(0 - < - Ha K ) + Ka K ){-±- - «)p(£)d£ 

Joo p^J 

Introducing the variable £ = l/£ we get the practical version 

r N = (N+ 1) £ K exp(iV + l)^^" 1 ) - «C _1 - $K) + «a«) (^py + P(CVC 

Results depend on k. This is the ordinary RG dependence that is strong enough 
because the accuracy is different at large momenta where we keep only the expansion and 
at small momenta where the exact formulae have to be used to make integrals finite. The 
obtained results are given in Table |]. Here a K is found from integration, a n s a%ve is found 
from the naive (to third order) relation between the schemes 

A-7T 

a « = _ K + Ka 3 ); 



while a e s xact is found from exact formulae for RG scheme relations (|T5| , [T6D . 

For finding the parameter a K from the integral it is useful to know the derivative of 
the integral with respect to the boundary value a K (m^) = a 

dr 1 

da p K (a ) 

We can find the derivative RG equation for coupling a K describing its dependence on 
the parameter k 

= a K p K (a K ). (27) 
This is a particular case of RG equations 

da n . . f a x n+3 dx 

that describe the dependence of the running coupling on coefficients of the (3 function. 
Note that the dependence on c is fixed by the choice of the parameter A to be A^g. The 
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0.368(16) 


0.380(18) 
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0.2209 
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0.311(07) 


0.319(08) 


1.9 


0.2037 


0.307(07) 


0.315(08) 


2.0 


0.2016 


0.304(07) 


0.313(08) 


2.1 


0.2003 


0.303 


0.312 


2.2 


0.1993 


0.303 


0.311 


2.3 


0.1986 


0.303 


0.312 


2.4 


0.1982 


0.303 


0.312 


2.5 


0.1979 


0.303 


0.313 


2.6 


0.1977 


0.304 


0.315 


2.7 


0.1975 


0.305 


0.316 


2.8 


0.1975 


0.306 


0.318 


2.9 


0.1975 


0.307 


0.319 


3.0 


0.1975 


0.308 


0.321 


3.2 


0.1976 


0.310 


0.325 


3.4 


0.1978 


0.313 


0.329 


3.6 


0.1980 


0.315 


0.334 



Table 1: k scheme results 
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K 




1.5 


.2058 


1.6 


.2049 


1.7 


.2039 


1.8 


.2030 


1.9 


.2021 


2.0 


.2012 


2.1 


.2003 


2.2 


.1994 


2.3 


.1986 


2.4 


.1978 


2.5 


.1969 


2.6 


.1961 


2.7 


.1953 



Table 2: RG dependence of a K on k 

extraction of a K is done under the assumption of RG invariance of the expression for the 
r lepton width so the only reliable data for a K can be taken from that part of Table 1 
where equation (^7]) is satisfied. This equation can be easily solved analytically (it is 
a linear equation if k is considered as a dependent function and a K as an independent 
variable), we have preferred however to solve it numerically in the vicinity of the value 
(k = 2.1, a K = 0.2003). The solution does not match well the pattern of the dependence of 
the extracted a K presented in Table 1 that means that higher order terms of perturbative 
expansion for the width are essential. If the theoretical expression for r T was exactly 
independent of the scheme the extracted charge would satisfy eq. (^). The actual solution 
with the initial condition (k, = 2.1,a K = 0.2003) which is our best estimate is presented 
in Table 2. So the correct pattern is around 1.9 < k < 2.3 where the prediction is 
a s (m 2 T ) = 0.312. 

Note that contrary to possible impression the prediction in fixed point scheme (or in K 
scheme) is also nonstable. Indeed, it is easy to introduce a set of schemes parameterized 
with the fixed point value of the invariant charge - the extracted a s (m 2 ) will depend on 
the scheme within the set. A {3 function for such a set could have the form 

0f(a) = P(a)(l + kP(o)) (28) 

that introduces a dependence on an external scheme parameter k. All schemes of the 
type fl2"8|) have a fixed point with different value of the coupling constant depending on 
parameter k. 



15 



To conclude, we have exploited the renormalization scheme freedom to show that some 
infinite subsets of diagrams can be even convergent not just summable in some generalized 
way. Then with lack of any parameter or criterion for choosing a particular set of diagrams 
the renormalization scheme freedom - extended to partially summed infinite series - can 
change the conclusion about the pattern of behavior of PT series in high orders dictated 
by the leading order RG and by the presence of the infrared pole in the running coupling 
constant. The very existence of renormalons can well be the attribute of a renormalization 
scheme with the Landau ghost. 

In fact, the dependence of the extracted values of the coupling constant for the r 
lepton width on schemes is rather large because the energy scale (m^) is quite low. So, 
it might be reasonable to fix the scheme in an arbitrary, and somehow simplest, way and 
then to parameterize the low momenta region in the integral sense only without detailed 
description of the behavior of the running coupling from point to point. This can be 
done in terms of distribution. Adding the localized distribution like 5 function and its 



derivatives can make the integral well defined |2(| These localized contributions look 
as nonperturbative terms. 



As for phenomenological applications, e.g. |32], [33], [34], [35|, p6| , there are models where 
some extrapolation of the running charge into the infrared region is used for practical 
calculation (normally under the sign of integration). In case of the r lepton width the 
uncertainties are large. Our results show that with an extrapolation chosen one can not 
freely add the contribution of gluon condensate to take into account nonperturbative 
effects and "improve" the computation - it must be coordinated with the continuation of 
the running coupling constant into infrared region. In fact, the fixed gluon condensate 
corresponds to a fixed continuation - the pattern of principal value looks very attractive 
in this respect. An analogous situation appears in a toy model considered in ref. [j3?|. 

Thus, the integration over an infrared region involves the strong coupling dynamics 
and is arbitrary to a large extent. The final expression of the results through the coupling 
constant at the boundary value is misleading - the series is bad defined in PT and requires 
some additional inputs concerning the behavior of the running coupling constant at small 
momenta. The operation of integration must be restricted, it remains valid only if the low 
s are suppressed by the weight function. Then the result of integration is perturbatively 
connected to the integrand and the procedure can be formulated in a RG invariant way. 

Acknowledgments. 

One of us (NVK) is indebted to participants of CERN Renormalon Workshop for very 
hot discussions that was one of the motivations for writing this paper. 



16 



References 



K.G.Chetyrkin, "New computing techniques in physics research III", 
eds. K.H. Becks and D.Perret-Gallix (World Scientific, 1993), p.523. 

G. 't Hooft, in: The Whys of Subnuclear Physics, Proceedings of the 15th Inter. School 
on Subnuclear Physics, Erice, Sicily, 1977, edited by A.Zichichi (Plenum Press, New 
York, 1979), p.943. 

F.David, Nucl.Phys. B209(1982)433. 

A.H.Mueller, Nucl.Phys. B250(1985)327. 

V.I.Zakharov, Nucl.Phys. B385(1992)452. 

M.Beneke, Phys.Lett. B307(1993)154; Nucl.Phys. B405(1993)424. 
P.Bali, M.Beneke and V.M. Braun, Nucl.Phys. B452(1995)563. 
M.Neubert, CERN-TH- 7524-94, Feb. 1995, pep-ph/9502264. 



D.J.Broadhurst and A.G.Grozin, Phys.Rev. D52(1995)4082; 
M.Beneke, V.Braun, Phys.Lett. B348(1995)513; 
C.N.Lovett-Turner, C.J.Maxwell, Nucl.Phys. B452(1995)188. 

R.Coquereaux, Phys.Rev. D23(1981)2276. 

G.Altarelli, PNason and G.Ridolfi, Z.Phys. C68(1995)257. 

A.L.Kataev and V.V.Starshenko, Mod.Phys.Lett. A10(1995)235. 

M.Neubert, CERN-TH-95-112, Oct. 1995, |hep-ph/9509432 . 



G.Grunberg, Phys.Lett. B221(1980)70; Phys.Rev. D29(1984)2315. 
N.V.Krasnikov, Nucl.Phys. B192(1981)497. 

A.L.Kataev, N.V.Krasnikov and A.A.Pivovarov, Phys.Lett. 107(1981)115; 
Nucl.Phys. B198(1982)508. 

N.V.Krasnikov, K.G.Chetyrkin and A.N.Tavkhelidze, Phys. Lett. 76B(1978)83; 
N.V.Krasnikov and A.A.Pivovarov, Phys.Lett. B112(1982)397; 
N.V.Krasnikov, A.A.Pivovarov and N.N.Tavkhelidze, Z.Phys. C19(1983)301. 

N.V.Krasnikov and A.A.Pivovarov, Phys.Lett. B116(1982)168. 

17 



A.A.Pivovarov, Z.Phys. C53(1992)461; Nuovo Cim. 105A(1992)813. 
N.V.Krasnikov and A.A.Pivovarov, |hep-ph/9510207 . 



O.V.Tarasov, A.A.Vladimirov, A.Yu.Zharkov, Phys.Lett. B93(1980)429. 

G.Parisi, Phys.Lett. 76B(1978)65; Nucl.Phys. B150(1979)163. 

E.Braaten, Phys.Rev.Lett. 60(1988)1606; 
S.Narison and A.Pich, Phys.Lett. B211(1988)183; 
E.Braaten, Phys.Rev. D39(1989)1458; 

E.Braaten, S.Narison and A.Pich, Nucl.Phys. B373(1992)581. 

S.G.Gorishny, A.L.Kataev and S.A.Larin, Phys.Lett. B259(1991)144. 

L.R.Surguladze, M.A.Samuel, Phys.Rev.Lett. 66(1991)560, 2416(erratum); 
Phys.Rev. D44(1991)1602. 

A.Dhar and V.Gupta, Phys.Rev. D29(1984)2822. 

S.J.Brodsky, G.P.Lepage and P.B.Mackenzie, Phys.Rev. D28(1983)228. 

M.Neubert, Phys.Rev. D51(1995)5924. 

P.M.Stevenson, Phys.Rev. D23(1981)2916. 

A.C.Mattingly and P.M.Stevenson, Phys.Rev. D49(1994)437. 

N.V.Krasnikov and A.A.Pivovarov, |hep-ph/9512213 . 



I.I.Bigi, M.A.Shifman, N.G.Uraltsev, A.I.Vainshtein, Phys.Rev. D50(1994)2234. 
M.Beneke and V.M.Braun, Nucl.Phys. B426(1994)301. 

A. V.Manohar and M.B.Wise, Phys.Lett. B344(1995)407. 

B. R.Webber, B339(1994)148; 

Yu.L.Dokshitzer and B.R.Webber, Phys.Lett. B352(1995)451. 



R.Akhoury and V.I.Zakharov, UM-TH-95-19, Jun 1995, |hep-ph/ 9507253 . 



A.A. Penin and A.A. Pivovarov, Preprint INR-894-95, Jun 1995. |hep-ph/950637C 



18 



